Spin Caloritronics in Noncondensed Bose Gases 



C. H. Wong,* H.J.van Driel, R. Kittinaradorn, H.T.C. Stoof and R.A. Duine 

Institute for Theoretical Physics, Utrecht University, 
Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 

We consider coupled spin and heat transport in a two-component, atomic Bose gas in the noncon- 
densed state. We find that the transport coefficients show a temperature dependence reflecting the 
bosonic enhancement of scattering, and discuss experimental signatures of the spin-heat coupling in 
spin accumulation and total dissipation. Inside the critical region of Bose-Einstein condensation, we 
find anomalous behavior of the transport coefficients, and in particular, an enhancement for the spin 
caloritronics figure of merit that determines the thermodynamic efficiency of spin-heat conversion. 



The Seebeck and Peltier effects are thermo-electric 
phenomena that are well understood for ordinary conduc- 
tors. Besides being theoretically interesting, these effects 
have many commercial applications ranging from wine 
coolers to thermo-electric generators. Recent develop- 
ments in spintronics have led to the nascent field of spin 
caloritronics [1] that introduces a spin-dependent gener- 
alization of these phenomena. In solid-state systems, the 
spin-Seebeck and spin-dependent Seebeck effects have re- 
cently been measured |2|l3j , while their theoretical expla- 
nations are still under active debate [4 . 

Cold-atom systems provide a perfectly clean environ- 
ment to study spin caloritronics, without many of the 
factors that can make the interpretation of solid-state 
experiments difficult. In addition, cold atoms are ad- 
vantageous in the study of spin and spin-resolved heat 
transport because one can experimentally realize systems 
where spin is conserved, apply different temperatures for 
different atomic species, and measure their distribution 
functions separately. While equilibrium properties of 
these systems have been thoroughly studied, recent re- 
search has started to focus on non-equilibrium behavior 
such as spin dynamics [5 , heat transport [6j, and spin 
drag |7|8j . The spin drag relaxation rate has been mea- 
sured in Fermi gases [8 , and experiments to measure 
the spin drag conductivity in Bose gases have recently 
been performed [9 However, the concomitant thermo- 
spin phenomena remain largely unexplored. 

In this Letter, we study the coupling of spin and heat 
transport in a cold atomic Bose mixture of two spin 
species and calculate the spin and heat transport co- 
efficients in the noncondensed state. In particular, we 
show that the bosonic nature of the particles leads to 
qualitatively different temperature dependence of these 
coefficients as compared to electronic systems. Further- 
more, we introduce a spin caloritronic figure of merit for 
this system called "Z^T", analogous to the "ZT" figure 
of merit that determines the efficiency of devices based 
on the usual solid-state thermo-electric effect. We com- 
puted the temperature dependence of Z^T and find an 
initial downturn on approach to the critical temperature 
of Bose-Einstein condensation, followed by an enhance- 
ment inside the critical region. This result is interesting, 
both in the context of developing high efficiency "atom- 
tronic" devices based on thermo-spin phenomena and the 
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FIG. 1: Schematic illustration of the coupled spin and heat 
transport. Spin-dependent forces or temperature gradients, 
labeled by F^,^., VT-^,^, generates both spin and spin heat cur- 
rents proportional to the coefficients Gs and Ps in former and 
ns and Ss in the latter case , as illustrated in the figure. 



study of dynamical critical phenomena [TO] Our theoret- 
ical results also have implications for thermo-spin phe- 
nomena in other systems where the transport is medi- 
ated by degenerate bosons, such as a quasi-equilibrium 
magnon gas 

We begin by considering a cold boson system above 
the critical temperature Tc of Bose-Einstein condensa- 
tion, composed of two different spin states selected from 
a larger integer- spin multiplet, which we will label "spin 
up" (t) and "spin down" {\). We apply forces and tem- 
perature gradients which are equal and opposite for the 
two spin species, i.e., = — F^ and VT-^ = — VT^ 
[12 . In response to the "spin force" and "spin temper- 
ature gradients" defined by Fg = F^ — F^ and VT^ = 
VT^ — VT^, there will be a spin current and spin heat 
current, js = jt - j;, and Qs = qt - Qi, respectively. We 
define the linear-response coefficients by (shown schemat- 
ically in Fig. [T]): 
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In the above, is the spin conductivity, k,'^ is the spin 
heat conductivity at zero F^, and and Pg are the 
spin-Seebeck and spin-Peltier coefficients, respectively, 
which are related by the Onsager reciprocity principle, 
Pg = SsT^ T being the temperature. Although we have 
written the response matrix in a form analogous to the 
thermo-electric coefficients defined in metals, we note 
that here the microscopic mechanisms for the coupling 
between spin and heat fiows are different than those in 
metals because there is no disorder and lattice phonons in 
the Bose gas. The spin- heat coupling studied here is akin 
to the thermo-diffusion effect in multicomponent classical 
gases [13 . The spin conductivity is, to the leading order, 
determined by the viscosity between up and down atoms 
that arises from inter-spin scattering, which is called spin 
drag. In contrast, the spin heat conductivity, which has 
dependence on intra-spin scattering, is finite even in the 
absence of inter-spin scattering. 

Using the Boltzmann equation for a two-component 
Bose gas, we have computed <Js, and as a function 
of nA^, with n the equilibrium particle density per spin 
state th at is assumed to be equal for both spin species and 
A = ^J2^T^l? /mkBT the thermal deBroglie wavelength, 
where m is the particle mass, h the Planck's constant, 
and Ub the Boltzmann constant. The results are shown 
in Fig. [2] (inset) and [sj The order of magnitude of the 
Seebeck coefficient is^^ ~ — ~ —1 /ieV/K, compa- 
rable to what is found in ferromagnetic materials [3] . No- 
table is the sharp decrease of all transport coefficients as 
one approaches the critical value nA^ 2.61. This effect 
is due to bosonic enhancement of scattering into occupied 
states, which dramatically increases as one approaches 
the critical temperature Tc = 27rh'^ /mkBA.^ where Ac ^ 
(2.61/n)^/^. This effect is illustrated in Fig. [i) in which 
we plot the effective relative momentum distribution /rei, 
defined so that the interspin collision rate is given by 

1/Tinter = CTinter / P^G^Pr^r/rel(Pr)/(27rfl)^, whcrC (Jinter = 

47ra^ is the interspin scattering cross section and Vr the 
relative velocity. As shown in the figure, this distribution 
increases sharply as one approaches the critical tempera- 
ture. This is in contrast to spin drag in degenerate Fermi 
gases, where due to Pauli blocking, below the Fermi tem- 
perature the spin conductivity increases as the tempera- 
ture is lowered [14 . 

Next we outline the calculation of the transport coeffi- 
cients outside the critical region. We start with the two- 
component, static Boltzmann equation for the semiclas- 
sical distribution function npcr(r) under spin-dependent 
external forces f^r. 



(vp • V + • Vp)np^ = Cp^[n^,n;], 



(2) 




FIG. 2: Spin diffusivity [normalized by (/i/m)(A^/a^)] as a 
function of the degeneracy parameter nA^. Inset: A plot 
of the spin conductivity as relative to the classical value 
(left axes) cts^^'* = (3 / 64^/27t) (A / a'^ h) and a plot of the spin 
heat conductivity at zero spin current, Ks, normalized by 
klTA/a'^h (right axis). 




where Vp = p/m is the particle velocity, a =t, i the 
(pseudo)spin of the two-component Bose gas, and the 



FIG. 3: Spin-Seebeck coefficient in units oi ks (right axis) 
and the dimensionless figure of merit ZsT (left axis). The 
plots here are outside the critical region, where the absolute 
magnitude of both coefficients decreases with temperature on 
the approach to the critical point (Note that the values of Ss 
are negative). 



collision integrals Ccr are given by 

[ni<^n2r(l + riza){l + n4r) - n3„n4r{l + ni„){l + n2r)]- 

(3) 

The collision integral describes the 2-body elastic scat- 
tering of particles labeled by (picr, p2r) {pso'^p^r), 
and is the solid angle between ingoing and outgoing 
relative momenta Pr = (pi— p2)/2 and p^ = (ps — P4)/2, 
respectively. We take the inter-spin differential cross sec- 
tion to be da^i/dQinter = ci'^, where a is the scattering 
length, and take the intra-spin terms to be da^^/dQ = 
da^i/dQ = 2a^ on account of Bose statistics ^6j. We 
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FIG. 4: A plot of an effective relative momentum occupation, 
/rei(Prei), for three values of the parameter nA^. The figure 
clearly shows a sharp increase in frei as one approaches the 
critical temperature. The Feynman diagram indicates that 
the interspin scattering rate is increased by the Bose enhance- 
ment factors, (1 + /-[)(! + /|), for the outgoing states. 



parametrize the non-equilibrium, steady state distribu- 
tion by 



.(r) = /p,(r)-(9e/>p-(r), 



(4) 



where = (exp[(ep — ii)/kBT] — is the equilib- 

rium Bose distribution, /i is the chemical potential, and 
fpai^.t) = (exp[ep - /i^(r))//cB^a(r)] - 1)""^ is the lo- 
cal equilibrium distribution, def^ = — /p(l + /pj/Zc^T, 
and (j)pa- describes the response to the spatial inhomo- 
geneities in /per and contain all dissipative effects. This 
parametrization represents an expansion in the ratio of 
the mean free path to spatial gradients [15]. 

Linearizing the Boltzmann equations with respect to 
(j)pcr and gradients in /pcr(r), we find that, for our choice 
of scattering lengths, the equation for the spin distribu- 
tion — n^i decouples from the equation for the total 
distribution -\- n^. We consider the linear response of 
the spin distribution to spin forces fg, gradients in the 
spin chemical potential /ig, and spin temperature Tg in- 
dependently from the response to the average forces, i.e., 
we consider the spin dependent forces = crfsf^ and 
VTa- = aVTs/2. The linearized drift terms in the left- 
hand side of Eq. ([2| are proportional to VTg and V/ig, 
but they are not independent thermodynamic forces be- 
cause the chemical potential has dependence on the den- 
sity and temperature, fia- = ficrina^T^j). We therefore 
transform the drift terms using the Gibbs-Duhem rela- 
tion, djiis = —sdTs + {l/n)dps-, and the thermodynamic 
identity w = jj. -\-Ts^ where s is the equilibrium entropy 
per particle, w is the enthalpy per particle, Ps = — Pi 
is the spin pressure, and we consider the response to the 
total thermodynamic force Fg = fg — ^Ps/n and VT^. 
Writing the linear-response solution as 



</.p, = <^F(p)-F, + (/.T(p)-(-Vr,), (5) 

the linearized collision integral in the spin equation can 
be written as Cps{(l)s) = Cpicf^r) ■ F, + Cp(0t) • (-VT,). 
The linearized Boltzmann equation for the spin distribu- 
tion requires that 4>f^ 4>t satisfy 

5J°Vp = Cp(0^), d^Jl (^^) Vp = Cp(0T). 

Using the solution to Eq. ([6|, the spin and spin heat 
currents are given by 



d^p 
d^p 



9efp Vp0p5 

dX (ep - ^)vp(/>p. (7) 



To solve Eq. ([g]), we expand the solutions in 
a power series, 0^ = Z)n=o ^^(m, ^) (ep/^B^)'' P, 
(pT = ^n=o^n{fJ^^T) {cp/kBT)^ and take moments 
of the Boltzmann equation by multiplying Eq. (|6| by 
(ep/Zc^Tj^p and integrating over p, resulting in a se- 
ries of equations of an^bn, which we truncate and solve 
at the second order. The transport coefficients are read- 
ily expressed in terms of the expansion coefficients and 
the temperature dependence were computed numerically. 
The Bose enhancement of the spin drag rate in the ab- 
sence of spin heat currents, calculated in Refs. |17|18j , 
was computed using the leading-order solution which de- 
scribes local Bose distributions for the spin up/down par- 
ticles rigidly shifted apart, resulting a spin current. The 
second-order solution which we have determined here 
represents a distortion of the local Bose distribution and 
is necessary to capture coupled spin and heat flows [19]. 

We note that (jg, which determines the spin current 
driven by external forces, is related to the spin diffu- 
sivity Dg^ which determines the spin current driven by 
diffusive forces induced by spin density gradients via 
is = —Dg^Us^ where = — is the spin den- 
sity. This diffusive current tends to return the system 
to homogeneous equilibrium and must satisfy the Ein- 
stein relation, Ds = (Js/Xsi where Xs = dns/djis^ is the 
static spin susceptibility and jUg = jii^ — f^i is the spin 
accumulation. We plot the spin diffusivity in Fig. |2] us- 
ing the noninteracting spin susceptibility. The diffusivity 
also determines the spin density gradient induced by a 
spin temperature gradient when the spin current is zero. 
For example, for a typical density of n = 10 
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cm 



at 



a temperature of T = 1.5 /iK (nA^ = 2.41), we find 
\\/ns/\/Ts\ = XsSs - 10^^ /ii^-^cm-^. For a typi- 
cal temperature gradient of 1 /iK/mm and for a cloud 
size of 1mm [6|, we find a spin density accumulation of 
Sug = 10^^ cm~^ , which gives a sizable experimental 
signal of a dus/n = 10%. 

Similarly, the heat current driven by VT^ tends to re- 
duce unequal distributions of energy between spin up and 
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down particles. The total dissipation is given by 

m5=^F2 + ^(VT,)2 + <T,5.F,-VT,, (8) 

where S is the total (spin-summed) entropy. It is thus 
possible to measure k.^ and Ss by measuring the heat- 
ing. Furthermore, the last term, analogous to Thomp- 
son heating, is sensitive to the relative signs of and 
VT5, and allows one to clearly distinguish experimentally 
the heating contribution from the spin-Seebeck coupling. 
It is important to note that typically, experiments are 
done in the presence of a trapping potential, which in- 
troduces spatial dependence in the transport coefficients. 
The measured values of the transport coefficients should 
be compared with the trap- averaged values, which dif- 
fers from the results presented here, but may readily be 
computed using our Boltzmann formalism. 

For ordinary conductors, one defines a dimensionless 
figure of merit ZT = cfS'^T/k^ with a the conductivity, 
n the heat conductivity in the absence of current and 
S the thermo-electric Seebeck coefficient, which deter- 
mines the efficiency of engines based on thermo-electric 
effects. Analogously, we define Z^T = cj^SIT / where 
— ^'s ~ ^sS^T is the spin heat conductivity at zero 
spin current, which is plotted in the inset of Fig. |2] [20 . 
We plot ZgT as a function of nA^ in Fig. [s] , and ob- 
serve an initial decrease as one approaches Bose-Einstein 
condensation. The quantity ZgT measures the ratio be- 
tween the magnitude of the spin heat coupling and total 
dissipation, thus the downturn observed results from the 
faster rate of decrease in the spin heat coupling compared 
to the total dissipation as one approaches Bose-Einstein 
condensation. 

However, the results based on the semiclassical Boltz- 
mann equation do not capture the critical fluctuations 
near the phase transition of Bose-Einstein condensation. 
It is known from the theory of dynamical critical phe- 
nomena that transport coefficients show anomalous be- 
havior in the critical region [10 which is experimen- 
tally accessible [21 . In this region, we have computed 
the transport coefficients using the Kubo formula with 



the Hamiltonian density 1-L = 



2m 



^^^^Trai^ti^ti^ri^ai whcrc ipcr diie the bosonic field 
operators and the strength of the contact interaction 
is given by the two-body T-matrix element: Terr = 



ATrh^acrr/'^' By neglecting the vertex corrections, the 
transport coefficients can be expressed in terms of the 
spectral function pcr(k, cj) = — Im[G^(k, cc;)]/7rfi, which 
is proportional to the imaginary part of the Fourier 
transform of the one-particle retarded Greens function, 
G^(r,t;r^tO = ie{t - t'){[ilj^{Y,t),^l)^ {y' ,t')]) . The spin 
conductivity then reads 122] 



l2m?kBT J (27r) 



sinh^(fiw/2fcBT) 



(9) 

while (JsSs and n'^ are given by including an additional 
factor (ep — ii)/T and (ep — /i)(ep — w)/T^ respectively, 
in the integrand above. Near the Bose-Einstein phase 
transition, according to Eq. ([9|, the power law scaling 
of (Js follows from the scaling equation for the spectral 
function, i.e., 

p(k, c^, /i - /i,) = A2-V(Ak, yuj, Ai/^(/i - /i,)) , (10) 



where A is an arbitrary scaling parameter and fic is the 
critical chemical potential. Furthermore, rj is the anoma- 
lous dimension, v is the exponent for the correlation 
length ^ - (/i - iic^ {T - TJ"^ and z is the dy- 
namical exponent. The theory of static and dynamical 
critical phenomenon [10 predicts that the exponents are 
close to their the mean field values of77 = 0, i^ = l/2, 
and z = 2, in agreement with our numerical calculations. 
With these values of the exponents and in three dimen- 
sions we find we find an enhancement of the figure of 
merit which scales as Z^T ~ ^ oc (T — Tc)~^/^. A numer- 
ical calculation [22 shows that the temperature range for 
this upturn is |T - Td/Tc ^ 60(a/Ac)^ Furthermore, in 
two dimensions, the critical region given by the Ginzburg 
criterion is larger, |T — Td/T^ ex a/ Ac, while the nature 
of the phase transition is different, and thus warrants fur- 
ther study. Lastly, since this result is generic to systems 
in the same universality class, we may expect similar be- 
havior in other systems of quantum degenerate bosons, 
such as quasi-equilibrium magnons. 
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